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Abstract 


Based  on  a  data  matrix  X  ■  (X^,.  .,X^):  pxk  with  independent  columns 

X  ~N  (5., I),  and  an  independent  Wishart  matrix  S:  p*p~W  (n,£),  testimators  dom- 
,1  p  _1  -  P 

inating  the  best  equi variant  estimators  of  jz  and  | £ |  are  obtained  under  two  types 
of  entropy  loss.  For  simultaneous  estimation  of  the  mean  vector  and  the  variance 
covariance  matrix  of  a  multinormal  population,  a  suitable  entropy  loss  is  developed 
and  testimators  dominating  the  pair  consisting  of  the  sample  mean  vector  and  the 
best  multiple  of  the  sample  Wishart  matrix  are  derived,  A  technique  of  SINHA 
(Jour,  Mult.  Analysis,  1976)  is  heavily  exploited. 


AMS  1980  Subject  Classification.  Primary  62F10;  Secondary  62H99 
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1.  INTRODUCTION.  Suppose  Y^,...,Yn  are  iid  N(£,cr2).  If  £  is  known,  then  the  best 

2 

scale  invariant  estimator  of  cr  is  given  by 

Cti-+*2)“1X"MlCYi-  O2.  (1.1) 

It  is  proved  in  Girshick  and  Savage  (1951)  ,  and  Hodges  and  Lehmann  (1951)  that 

2 

is  an  admissible  estimator  of  cr  under  squared  error  loss.  However,  if  £  is  unknown, 
then  Stein  (1964)  has  shown  that  the  natural  estimator 

<|.1(Y1,...,Yn)-  (n+l)'1^_1CYi- Y)2,  (1.2) 

of  a2(Y*  n  ^FJ-jY^)  is  inadmissible  under  squared  error  loss,  and  is  dominated  by 
estimators  of  the  form 

HY1,...,Yn)-min((n+l)'1^=1(Yi-Y)2,(n+2)“1^>tl(Y1-50)2]  (1.3) 

2 

for  every  fixed  constant  5q.  The  estimator  <J>  of  a  can  be  viewed  as  a  preliminary 

“In  “2  2 

test  estimator  (testimator)  which  uses  the  estimator  (n+1)  £^=^(Y^-Y)  for  cr 

if  the  F-statistic  n(Y-  Cq) { (n— 1)  ^^_^(Y^-Y) 2 }  for  testing  Hq :  5 =  5q  against 
the  alternatives  H^:  £  t  5^  exceeds  (n-l)/(n+l)  (thereby  rejecting  H^  at  a  certain 
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*  1  n  2 

significance  level),  and  uses  the  estimator  (n+2)~  ^.^(Y^Cq)  otherwise.  Stein 
(1964)  considered  the  more  general  regression  analog  of  this  problem  in  the  canonical 
set  up.  Brewster  and  Zidek  (1974)  have  shown  that  the  results  extend  to  a  more 
general  loss  including  the  entropy  loss  (first  introduced  in  James  and  Stein  (1961)) 
given  by 

L(a,cr2)  =  a/a2  -  log(a/a2)  -  1,  (1.4) 

to  which  attention  will  be  restricted  in  this  paper. 

There  are  two  possible  multivariate  extensions  of  the  above  results.  One  can 

consider  estimation  of  the  variance-covariance  matrix!  or  the  generalized  variance 

|e|  in  a  multinormal  set  up.  To  fix  ideas,  let  Y, ,...,Y  be  iid  N(£,E),  where  each 
~  -1  _m  „  _ 

is  p*l.  When  both  £  and  E  are  unknown,  the  minimal  sufficient  statistic  for  these 

parameters  is  (X,S) ,  where  X  -  m1/2Y  and  S  -  7“  , (Y  -Y  )  (Y  -Y  )T(Y  * m-1y“  ,Y4). 

~  ~  -  _m  „  i=I  _l  _m  _i  _m  °i=l_i 

Haff  (1979b,  1980,  1982)  and  Dey  and  Srinivasan  (1985)  have  considered  estimation 
of  E  under  several  losses  including  the  entropy  loss 

LjCE.E)  -  tr(EE_1)  -  logjEE"1!  -  p.  (1.5) 

They  propose  estimators  E  of  E  which  are  functions  of  the  Wishart  matrix  S  alone, 
but  do  not  consider  any  Stein-type  estimators  (i.e.  testimators) . 

In  this  note,  we  consider  estimation  of  E  and  E  ^  each  under  the  entropy  losses 
and 

L2(E,E)  =  tr(EE-1)  -  log| EE_1|  -  p.  (1.6) 

To  our  knowledge,  the  loss  (1.6)  has  never  been  considered  before  either  for 
estimating  E  or  I  \  Haff  (1977,  1979a,  1979b)  considers  estimation  of  E  *  under 
various  quadratic  losses.  For  us,  the  loss  (1.6)  seems  to  be  as  natural  as  (1.5), 
and  can  be  motivated  as  follows.  Suppose  S  is  a  random  variable  having  a  p- 
dimensional  Wishart  distribution  with  degrees  of  freedom  n  and  parameter  E(to  be 
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denoted  by  Wp(n,I)).  Write  fj.(s)  as  the  pdf  of  S.  Then,  a  meaningful  loss  in  es¬ 
timating  Z  by  A  (or  Z  ^  by  A  is  the  entropy  distance  between  W^Cn.Z)  and  W^(n,A), 


and  is  given  by 


Ez[log 


fz(s) 

YXs T]  =  (n/2)L2(A,Z). 
A  _ 


(1.7) 


Use  of  an  estimator  Z  in  place  of  A  gives  rise  to  (1.6). 

In  Section  2  of  this  note,  we  consider  estimation  of  Z  and  Z  ^  each  under  the 
losses  (1.5)  and  (1.6),  and  develop  Stein-type  testimators  dominating  the  best 
multiples  of  the  Wishart  matrix  and  its  inverse.  A  technique  of  Sinha  (1976)  is 
heavily  exploited.  Incidentally,  it  may  be  remarked  that  for  the  loss  (1.6)  no 
Haff -type  improved  estimator  over  the  best  equivariant  estimator  is  readily  available. 
We  also  consider  simultaneous  estimation  of  the  mean  vector  and  the  variance- 
covariance  matrix,  and  develop  certain  estimators  dominating  the  pair  consisting 
of  the  sample  mean  vector  and  the  best  multiple  of  the  Wishart  matrix  under  a 
suitable  entropy  loss  to  be  developed  in  Section  2. 

In  Section  3,  we  consider  estimation  of  (z|.  This  problem  has  received  attention 
in  Shorrock  and  Zidek  (1976),  and  Sinha  (1976).  In  these  papers  Stein-type  testima¬ 
tors  are  developed,  and  are  shown  to  dominate  the  best  multiple  of  |s|  under  squared 
error  loss.  Similar  testimators  are  developed  in  Section  3,  and  are  shown  to  dominate 
the  best  multiple  of  |sj  under  the  two  entropy  losses 


Lx( | £ | , | Z  j )  *  | Z  | / | Z  |  -  log(|z|/|z|)-  1  (1.8) 

and 

L2 ( |  £  |  ,  |  Z  l )  -  |  Z  |  /  |  £  |  -  log(  |  Z  |  /  |  £  | )  -  1 .  (1.9) 

Throughout  this  paper,  for  two  matrices  A  and  B  of  the  same  order,  A>_B  implies 
that  A-  B  is  nonnegative  definite.  In  the  remainder  of  this  section,  we  state  with- 
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out  proof  three  matrix  lemmas  which  are  used  repeatedly  in  Section  2.  The  proofs 
of  these  lemmas  are  quite  straightforward. 

LEMMA  1.  Let  denote  the  class  of  all  nonsingular  matrices.  Then  for  any 

A  €  F  and  Be  F  , 

P  P 

tr(AB)  -  log  I B  j  >_  log  j  Aj  +  p,  (1.10) 

equality  holding  iff  B  *»  A 

LEMMA  2.  Suppose  A^  0  and  B  C,  where  A,  B,  C  and  the  null  matrix  0  are  square 
matrices  of  the  same  order.  Then, 

tr  AB>_  tr  AC.  (1.11) 

LEMMA  3.  For  any  positive  definite  matrix  A, 

tr  A  -  log|A|  -  p^.0, 

with  equality  iff  A  ■  I  . 


2.  ESTIMATION  OF  E  AND  E  \  Consider  a  multivariate  normal  linear  model  in  its 
canonical  form.  Suppose  X  *»  (X^,...,X^)  is  a  pxk  matrix  with  independent  columns 
Xf~  Np(^i,E) ,  and  let  S  be  a  p-dimensional  Wishart  matrix  with  degrees  of  freedom  n 
and  parameter  E  distributed  independently  of  X.  We  assume  n>  p+1  and  unknown. 

Consider  first  estimation  of  E  under  the  loss  (1.5) .  As  pointed  out  by  Shorrock 
and  Zidek  (1976)  ,  the  above  problem  remains  invariant  under  the  full  affine  group 
G  acting  on  the  space  of  pxk  matrices  (writing  E,  =  (£^» , . .  ,C^) ) 


X  -*-AX  +  B,  Z+  AS  +  B, 


S  -*■  ASA^, 


E  -*■  AEAT, 


(2.1) 
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where  A  is  any  nonsingular  p*p  matrix,  and  B  is  any  p*k  matrix.  Then,  any  affine 
equivariant  estimator  of  E  must  be  of  the  form 


<j>0(S)  =  cS , 


(2.2) 


where  c  is  a  constant.  Noting  that  E(S)  =  n  ,  it  follows  that  under  the  loss  (1.5), 

the  optimal  choice  of  c  minimizing  the  risk  of  cS  under  the  loss  (1.5)  is  c»n 

T  -1 

Following  Sinha  (1976) ,  write  S  =  WW  and  U  =  W  X,  where  W  is  a  p*p  non¬ 
singular  matrix.  In  order  to  improve  on  the  best  affine  equivariant  estimator 

-1  *  t 

n  S,  consider  the  class  C  of  estimators  of  E  having  the  form  E  =  <1>(W,U)  =  , 

where  ip  =  i|/(UU  )  is  a  pxp  nonsingular  matrix.  This  class  C  contains  estimators 

equivariant  under  a  nonnormal  subgroup  H  of  G  obtained  from  G  by  putting  B*  0. 

-1  T 

The  special  choice  \p  =  ipn  ■  n  I  leads  to  the  corresponding  estimator  ■  Wi|>nW  = 

«,  ,0  «p  ~  u  ~  u~ 


n  of  E.  In  order  to  compute  the  risk  R,  of  $  under  the  loss  (1.5),  first  let 

t  -  T 

X  *  AX,  W  *  AW,  where  A  is  a  p*p  nonsingular  matrix  such  that  AEA  -  I,  1  de- 

noting  the  identity  matrix  of  order  p.  Note  that  U*  W  X  =  X^.  Then,  writing 

E  as  the  expectation  under  N  (5. ,1)  for  X  i^  l,,,,,p  and  W  (n,E)  for  S, 

£  t  L  P  •  -**1  ** 

and"  5*=A5,  one  gets 

•*  -»n  _  . 

K  -  Er  -[trOtywV1)-  log|W^WTE"  |  -  p] 

<P  Z,1  - -  - - 


E  z  [tr(wjwj,  )  -  log|W*Wj  -  log  I*  |  -  p]. 


(2.3) 


Note  that  for  comparing  the  risk  performance  of  members  within  the  class 
Ca>  -  W*WT  (under  5,2)  *  (under  5^,1^, it  suffices  to  consider 


R,  »  E.  T  [tr(W*W  ^  )-  log S ^  1  1 


E[  tr{i|»~1(C*»U)  -  log  |  |  ] , 


(2.4) 
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where  «  E  (W^W*|  U=  u),  and  E  denotes  expectation  over  the  marginal 

'  ~  S  JL  I  i  ~  ~  ~  - 

~  ~P 

distribution  of  U. 

To  minimize  R  with  respect  to  <p ,  it  suffices  to  minimize 


tr{'l>*1(£*,u)'J<}  -  log  | 


(2.5) 


with  respect  to  ip  for  every  u.  Using  Lemma  1  with  A=  ip*  (5*>u)  and  B=  ip,  it  follows 
that  the  expression  in  (2.5)  is  minimized  when  tp  -  ip^CZ*, u) .  However,  this  expression 
involves  not  only  u  but  also.  We  find  next  an  upper  bound  for  i^(£*»u)  free 
from  £;*• 

LEMMA  4.  K(£  ,u)  <  t|>(u)  -  (n+k)_1(I  +uuT)  . 

-  ,  —  ..  _ p _ 

The  proof  of  this  lemma  is  deferred  to  the  Appendix.  Based  on  iKu) ,  a  testimator 

is  now  constructed  as  follows. 

~  t  _l  T 

Let  4>=WiJ/W  **  (n+k)  (S+XX  ).  For  estimating  E,  define  the  testimator 


4>  *  <(>  if  4  £  <PQ 


otherwise. 


(2.6) 


The  corresponding  ip  say  ip  is  given  by 


ip  ■  *P  if  P  1  'Pg 


ip^  otherwise. 


(2.7) 


Remark  1.  The  estimator  <i  defined  in  (2.6)  is  a  multivariate  generalization  of 
Stein's  (1964)  univariate  testimator.  The  condition  4><_4>q  can  be  alternately  ex¬ 
pressed  as  XXT<_  (k/n)S  <=>  sup(  2.TXXT£) <_  k/n  <->  largest  eigenvalue  of 

. 

T  —1  **  “ 

X  S  X<  k/n.  Thus,  the  estimator  proposed  in  (2.6)  is  based  on  Roy's  maximum  root 
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2 

test.  The  test  reduces  for  k-  1  to  Hotelling's  T  test.  I 

~  T  ! 

Remark  2.  The  condition  i|»  <_  can  be  alternately  expressed  as  UU  <_  k/n  <=>  ||  u||  <  k/nj 

where  |j  *||  denotes  the  Euclidean  norm. 

The  following  theorem  shows  the  dominance  of  the  testimator  p  over  <f>g. 

THEOREM  1.  Under  the  loss  (1.5), 

R~  <  R  for  all  £  and  E. 

-  ~ 

Proof:  Using  (2.3),  (2.4)  and  (2.7), 


where  in  the  last  step  of  (2.9),  one  uses  Lemma  3.  The  proof  of  Theorem  1  is 
complete  from  (2.8)  and  (2.9). 


Remark  3.  Quite  generally,  given  any  estimator  cp  *  WipW  ,  defining  p  =  4>  if  <p  <_  p , 

~  T 

*  otherwise,  where  <{i  *  W>pU  and  ip  is  defined  in  Lemma  4,  one  gets  R„  <  R  for  all 


5  and  E  (vide  Sinha  (1976)).  This  enables  one  to  develop  sequential  testimators 


as  in  Sinha  (1976). 
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Next  we  consider  the  loss  given  in  (1.6).  In  this  case,  the  best  multiple 

of  S  (minimizing  the  risk)  is  given  by  (n-p-1)  \  Let  ^qq(S)  =  (n-p-1)  , 

T  -1 

=  WiJ/.-W  so  that  _  =  (n-p-1)  I  .  Once  again,  we  consider  a  competing  class 
~  UU^  ^  UU  ^  p 

T  T 

Cq  of  estimators  of  the  form  <j>(S)  =  W^(UU  )W  .  Proceeding  as  in  (2.3),  under  the 
loss  (1.6)  ,  the  risk  of  4>  is  given  by 

R  =  Er  [tr((W*W*  fVVlogj^W*  )_1|  -  log  |*_1|  -  p]. 

T  ~  -  - 

~  ~P 


Hence,  for  comparing  estimators  of  the  given  type  <p,  it  suffices  to  consider 


R  »E  [tr{(WTW.  fV1}-  log|^_1|]. 

t  5*^P 


Using  Lemma  1  once  again,  it  follows  that  the  optimal  choice  of  ij;  is  E 

5** 

=  '^(-*tu)  (say).  Similar  to  Lemma  4,  we  now  prove  the  following  lemma. 


I 

-P 


(2.10) 

uwIVVi 


LEMMA  5.  i|»1(^J.,u)<  (n-p-l+k)  ^(I  +uuT)  =  ^(u)  (say). 

-  p  „  _ 

-  t 

The  proof  of  Lemma  5  is  also  deferred  to  the  appendix.  Let  <J>q (S)  =  W^^(u) W  . 
Similar  to  the  previous  situation,  we  define  the  testimator 

*0(s)  =  d»0(s)  if  d>0(s)  <_ 4>00 (s) 

=  tpQQ  (S )  otherwise.  (2.11) 


Accordingly,  the  ^  value  corresponding  to  4>q  will  be  given  by 

||q(u)  =  ^q(u)  if  ^q(u)  <_  J00(u) 

*>  'i'QQ(u)  otherwise.  (2.12) 

~  -IT 

Remark  4.  Note  that  4> q ( S )  =  (n-p-l+k)  (S+XX  ).  Hence,  the  condition  $q(S)  <_  ^^(S) 
can  be  equivalently  expressed  as  largest  eigenvalue  of  X'S  ^X^_ k/ (n-p-1) .  Thus, 
in  this  case  also,  the  preliminary  test  is  based  on  Roy's  maximum  root  test. 
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We  now  prove  the  following  theorem. 


THEOREM  2.  Under  the  loss  (1.6),  R_  <  R 


!0 


roo 


Proof:  Write 


R~  -  R  -  R.  -  K 

r  4>  7  ^ 

to  ~00  *0 


00 


I*-1! 


■EUrt^CC.-UX*'1-*^)}-  log  -r  J‘[?0<*00] 


]I 


(2.13) 


‘roo1 


In  view  of  Lemma  5,  putting  A  =  i p  ^  -<p^,  B  =  and  C  =  \p  (E.,U)  in  Lemma  2  one 

^  -.U  .*00  ~  -U  ^  ^1  .**  .* 


gets 


rhs  of  (2.13) 

,~-l  -1 


<  Eltr^Oj;’1-  tc0i)+i°g|t0t00|;iIti  <^nn] 


=  -E[tr(^QipQQ)  -  logltotoo I  ”  P]If  ~ 
<  0, 


iT  -li 


o-:ooJ 


[to^too] 


(2.1^) 


where  in  the  last  step  of  (2.14),  one  uses  Lemma  3.  The  proof  of  the  theorem  is 
complete  from  (2.13)  and  (2.14). 

Remark  5.  Here  again,  as  explained  in  Remark  3,  one  can  develop  sequential  testima- 


tors  each  dominating  the  best  equivariant  estimator  ^^^(S). 


We  consider  now  the  simultaneous  estimation  of  the  mean  vector  and  the  variance- 


covariance  matrix  under  entropy  loss.  Writing  f^  ^(Y^, . . . ,Y^)  as  the  joint  pdf  of  n 


iid  N  (u,E)  variables,  it  follows  that  taking  the  loss  L((u,E),  (l, A))  as  the  entropy 
P  -  ~  -  -  ~  - 


distance  between  the  N^(u,E)  and  Np(£,A)  distributions,  and  assuming  n>  p+1  for 


purposes  of  estimation,  one  gets 
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f  _(Y.  Y  ) 

y,Z  .1  _n 

L[(y,Z),(£,A)]  =  E  [log  ~  — 

-  ~  '  iJ»“  r£,A  .1 . ,n; 


n 

E  [log(  1 S  1  /  |  A l )  ^trCZ^S  )  +  n(Y-y)TZ_1(Y-y)} 


+  ji  tr  (A-1S0)  +  n(Y  -  l)  TZ_1(Y  -  £)  }  ] 

(where  Y  =  SQ  =  ^(Y^Y)  (Y±-Y) T) 

=  log|EA  1|  2  “  ^  "  2 

+  |tr(ZA_1)  +  i[p  +  n(H-y)Tr1C£-y)] 

I 

=  f[U-y  ?E_1(£-y  )+  tr(ZA_1  )-  log|ZA_1|-  p].  (2.15) 

Thus  if  one  uses  the  best  location  and  scale  invariant  estimators  Y  and  (n-p-1)  XSq 
for  y  and  Z  respectively,  in  view  of  the  loss  (2.15),  it  suffices  to  improve  on  Y 
and  (n-p-1)  separately.  From  James  and  Stein  (1961),  it  follows  that  under  the 
loss 


LQU,y)  -  (£-y)Xz'1U-y), 

Y  is  improved  by  the  estimator  <5n(Y,S-.)  =  (1 - _)Y  ,  provided  p>  2. 

-U  ~  ~U  nfsj  Y 

Again  the  loss 

L2(Z,A)  =  tr(ZA-1)-  log | ZA_1 |  -  p 

is  the  same  as  the  loss  (1.6),  and  taking  k*  1  in  Theorem  2,  it  follows  that 
(n-p-1)  XSq  is  dominated  by  the  testimator 

6X(Y,S0)  -  (n-p)"1(S0+nYYT)  if  (n-p) ^(Sq-KiYY1)  <  (n-p-1) ‘1S0 

-  (n-p-l)"1^ 


otherwise. 
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Thus,  under  the  loss  (2.15),  (Y,Sq)  is  dominated  by  (6g,6^)  for  p>  2.  For  p^2, 

(Y,Sq)  is  of  course  dominated  by  (Y,6^).  Finally,  for  p^_2,  sequential  testimators 

each  dominating  (Y,Sq)  for  the  loss  (2.15)  can  be  easily  obtained  ( vide  Sinha  (1976)). 

Next  we  consider  estimation  of  the  precision  matrix  I  ^  under  the  losses  (1.5) 

and  (1.6)  (calling  L^(E,E)  and  L^CZ.Z)  as  L^(Z  \z  ^)  and  I^E  \z”^)  respectively). 

T  -1  -1 

Consider  the  class  of  estimators  of  the  form  4>  =  (WtpW  )  for  Z  so  that,  under  L^, 

the  choice  ib  *  n  ^1  leads  to  the  best  multiple  (of  S  ^)  estimator  nS  ^  of  Z  ^ . 

~P  -  - 

"-IT  ~ 

Recalling  that  ( S)  *  (n+k)  (S+XX  ),  it  follows  that  defining  (j>  as  in  (2.6),  the 

best  equivariant  estimator  nS  ^  =  <j>Q^(S)  is  dominated  by  $  ^(S)  .  Similarly,  under 

the  loss  L^,  the  best  equivariant  estimator  (n-p-l)S  ^  *  4>qq(S)  of  Z  ^  is  dominated 

by  ^q^(S)  defined  in  (2.11).  Again,  in  each  case  sequential  testimators  are  easily 

obtained. 

3 .  ESTIMATION  OF  |z|.  Consider  the  same  set  up  as  of  Section  2.  Estimating  [z| 
by  a,  assume  the  loss  to  be  given  by 

L.  (a,  | Z  | )  -  — - log  — - 1.  (3.1) 

1  '  M  |i| 

Following  Shorrock  and  Zidek  (1976)  and  Sinha  (1976),  it  follows  that  the  best  equi¬ 
variant  estimator  of  |z|  is  Cq|s|  where  Cq  is  determined  from  minimizing  F.^  ^ 

-  ~P 

( c | S |  -  logc)  with  respect  to  c.  This  gives 

co  “  ^Ez-  i  ls^_1  “  (n-p)'/n*  • 

-  -P 

Following  Stein’s  suggestion,  and  arguments  as  in  Shorrock  and  Zidek  (1976) 
or  Sinha  (1976),  we  look  for  better  estimators  in  the  class  $(X,S)*  MX^S  ^X)|s| 


for  some  real  valued  function  \p.  Under  the  loss  (3.1),  <p  has  the  risk 


K  *  Er  t  t^(XTS_1X)  is|  -  log^(XTs"1X)-  log  |  S |  —  1],  (3.2) 

-  -P 

T  -1 

where  is  defined  as  in  Section  2.  Write  V  *  X  S  X  so  that  given  V  =  v,  the 
best  choice  of  \p  (minimizing  (3.2))  is  given  by  ip(v)  =  ip  (v)  *  (E  _  (|s|  V=v)}  \ 

5*  ~  ~**Zp  '  ~  " 

Following  the  line  of  argument  of  Sinha  (1976) ,  one  can  easily  show  that 

^  (v)i  |  I^+v)  (n-p+k)  !/(n+k) !  =  4>0(v)  (say).  (3.3) 

Then,  it  is  easy  to  show  using  (3.3)  and  strict  convexity  of  the  loss  (3.1)  that 
for  every  defining  i|>(v)  =  min(i}i(v)  ,4»q(v)  ) ,  i|>(X^S  ^X)|s|  dominates  ij>(XTS  ^X)|s| 
under  the  loss  (3.1).  In  particular,  the  estimator  Z »  min{  -|  S |  , |  S+XX^ |  } 
dominates  { (n-p) ! /n! } |s |  under  the  loss  (3.1).  Note  that  Z  is  indeed  a  testimator 
since  the  ratio  |s+XXT|/|s|  ■  |  I^+XTS  ^x|  is  a  MANOVA  test  statistic  for  H^:  5=0 
against  H^:  5^0, 

For  the  other  loss  l^a,  |  E  | )  defined  by 

,  ,  M  1*1 

L2(a,  | Z | )  =  — - log  - - 1  (3.4) 

it  follows  that  the  best  equivariant  estimator  of  |l|is  c|s|  where  c  minimizes 
K  i  C-^TgT  +  logc) .  This  gives  c  -  EJ,=1  {|s[  1}  =  (n-p-2) ! /(n-2) ! 

T  "1 

As  before,  we  look  for  a  better  estimator  in  the  class  <J>(X,S)  ■  <HX  S  X)  |  S  | 
for  some  real  valued  function  <Ji.  Suhh  a  <fr,  under  the  loss  (3.4),  has  the  risk 

R  *  E  _  (1/(^(XTS"1X)  |s|)  +  logiKXTS_1X)  +  log | S |  -  1]  (3.5) 

which  is  minimized  for  a  given  V- v  by  choosing  ij>(v)  *  ’!>-  (v)  =  E  T  {|s|  1|v«v}. 

'  '  -  5*  ’  p  **  "  ~ 

Following  Sinha  (1976),  we  can  easily  show  that 


i t>e  (v)  <  I.+v|  (n-p-2+k)  !  /  (n-2+k)  !  -  <Pn(v)  (say). 


(3.6) 
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Then,  for  every  <J>  defining  <Kv)  -  min(y>(v) ,  <pQ(v)),  it  follows  that  iKXTs"1X)jsj 

T  —1  "** 

dominates  tp (X  S  X)  |  S  J  under  the  loss  (3.4).  In  particular,  the  testimator 

Z  -  min{-gl£npl  |  s  |  ,  ^EZ^IL|s^xxT|}  dominates  {(n-p-2) !  / ( n— 2)  |s|  under  L2 

loss. 

As  in  Sinha  (1976),  it  is  possible  to  easily  derive  sequential  testimators  of  |  E { 
under  both  the  losses.  Details  are  omitted. 


APPENDIX 


Proof  of  Lemma  4.  Since,  I  +  uu  is  positive  definite,  there  exists  a  nonsingular  Q 

T  T  ~ 

such  that  QQ  ■  I  +  uu  .  Write  *  W^Q.  Then  the  inequality  i^A(CA,u)  <_  <Ku)  <=> 

~P  ~  ~  ~  _i  -  _i 

—  I  **-]_  T  T  T  t 

(u)  can  be  alternately  expressed  as  Q  E  T  (w**w**|u)Q  ^  (n+k)Q  Q 
~  ~  ~  •»  -  -  -  ~  ~  ~  - 

-  -P 

Hence,  it  suffices  to  show  that 


E5* ,Ip< '**-** (n+lt)~p  for  a11 


(A.  1) 


Note  that  (A.l)  can  be  alternately  expressed  as 


E  I  (£  w**w**£|u)  1  (n+k)lTl 
~*’~P  . 


(A.  2) 


for  all  Z(t  0),  u  and  From  (2.19)  of  Sinha  (1976),  it  follows  that  a  sufficient 


condition  for  (A. 2)  to  hold  is  that 

n^+k 


2  IWW1 1  2  (£TvTwi)exp£-  -|<wwT  -  2wu^S^)]dw 


n-p+k 

p2  |wwT|  2  exp[-  y(wwT  -  2wu^£^) ]dw 


>(n+k)it  Z 


(A.  3) 


for  all  Z^O,  u^Q  u  and  In  (A. 3)  and  in  what  follows,  we  use  the  notation 

w  for  wAJk  (and  accordingly  W  for  W^A). 

T  T 

Use  now  the  transformation  Z  =  WL  ,  where  L  is  an  orthogonal  matrix  with  its 

T  1/2 

first  column  vector  equal  to  Z/(Z  Z)  .  We  write  Z  ■  (Z^,...,Z  ).  Then  (A. 3) 
can  be  alternately  expressed  as 
,  n~Ptk 

I  2IzzT|  2  (z^z^expl-  7tr(zzT"2zuL5*)  ldz 


n-p+k 
T  i  2 


IT  T 

exp[-  -rtr(zz  -  2zuL£^)]dz 


>  n+k,  where  u.  *  Lu.. 


(A. 4) 


Next,  as  Sinha  (1976),  let  A^  ■  1^  and 

ZJ*1Z(4.1  •  •  •  ZLlZ  ]  fz,., 

*t  - 1- z„)  -1+l-1+1  -i+l-p  :1+1 

z  z. ..  .  .  .  zLz  ; 

'-.p.i+1  -p.p  '  |^Zp 

Is  !,•••, p—  1  • 


Then,  following  (2,20)-(2. 21)  of  Sinha  (1976)  and  noting  that  his  w's  are  our  z's, 
we  can  express  the  left  hand  side  of  (A.4),  in  Sinha's  (1976)  notation,  as 


n-p+k 


BlE* 

.  2  ,..T, 


n-p+k 
.  2  T 


E[(“pV  2  (ViWh’  •  •  ' 


n-p+k 


n-p+k 


n-p+k 


r  E[(W^W  )  2  (W*  A  W  )  .2..(«V«  )  2  ] 

-P-P  „p-l-p-l.p-l  ,1-1.1 


(A.  5  ) 


where,  given  , .  . .  ,W^  ^  ,  W  A^W^  is  a  noncentral  chisquared  variable  with  i.d.f. 

2  T 

and  noncentrality  parameter  -  n(i)Ain(i)  where  (i)  » *  *  *  *  n(p)  ^  ”  ul5** 

T  T 

Accordingly,  using  the  fact  that  >  W^A^W^  and  proceeding  as  in  Sinha  (1976),  we 
get  the  expression  in  (A. 5)  >_  2  (°  ^ffi)  +  p  -  n+k,  where  in  the  ultimate  step,  one 
uses  (2.22)  of  Sinha  with  r  ■  (n-p+k) /2.  The  proof  of  Lemma  4  is  complete. 


Proof  of  Lemma  5.  It  suffices  to  show  that 


E  [ lT  (W*W^)_1 1  |  u ]  <  ( n-p-  1+k)  -1  £T  ( Ip+uuT)  £ 


(A. 6) 


for  all  £(/0) ,  u  and  Defining  Q  as  in  the  proof  of  Lemma  4,  and  using  calculations 

similar  to  (A.1)-(A.3),  we  find  that  (A. 6)  can  be  equivalently  expressed  as 
n-p+k 

2|ww^|  2  (S,T(w^w)  ^i)exp[-  jtr(ww^-2wuAC^)dw 

wc  EP  , 

n-lH-k 

*  2)wwT|  2  expt-  ytr(wwT-2wu*S*) ]dw 


<  (n-p-l+k)  ^(£T£) 


(A.  7) 


"t  ^  ^**^7  v*  *  • » *7*  w » v  ■ 
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1  T 

for  all  l(.t  0),  u  -  Q  u  and  5..  Next  make  the  transformation  Z  *  WL  where 

*  ^  «  ■<  **» 

Z  -  (Z  ,..,z  )  and  L  is  an  orthogonal  matrix  with  its  first  column  vector  equal  to 

—  «*1  —  P 

H/(iT l)l,2m  Then,  (ZZT)-1  *  (lVwL)"1  =  LT(WTW)_1L  so  that  aT(WTW)-1y  uV)  is 

T  —1  T  -“1 

the  element  in  the  first  row  and  first  column  of  (ZZ  )  .  We  denote  this  by  (ZZ 

T 

Now,  writing  ^  -  L  u*,  (A. 7)  can  be  equivalently  expressed  as 
n-p+k 


[  |zzT|  2  (zz1)"1  exp[-  itr(zzT-2zTu^5*)]dz 

U-  ''  -  ’  . 


n-p+k 

|zzT|  2  exp[-  kzzT-2zTuT5*)]dz 
2  — 


<_  (n-p-l+k)  . 


(A.  8) 


With  the  same  A. ’s  (i=  l,...,p)  as  defined  in  Lemma  4,  it  follows  that 
*»  1 

I  I  -  ^p-lfp-lfp-1^  *  *  *  (?2^2?2)  ^l^lfl? 


and 


-  i??TrVpZp>  • 


•  (?2*2?2)- 


Accordingly,  writing  r  =  (n-p+k)/2,  lhs  of  (A. 8) 

.T .  .  .  r  f7^'A  7  >r  rn^x  •»  Nr~in 
^pAp)  ...  ^2a2z2; 

.T.  „  \r  , 


*  E[(?p-p?p)  •**  (?2^2~2^  ^l-l-l^  1 


-P-P-P 

<  (2r-l)_1(using  conditional  argument  as  in  the  proof  of  Lemma  4  and  (2.22) 
of  Sinha  (1976)) 

=  (n-p+k- 1)  ■*". 


The  proof  of  Lemma  5  is  complete. 
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